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SUMMARY 

Charts are presented that enable convenient determination of the 
pressure drop sustained by monatomic gases (ratio of specific heats, 
5/3) flowing at high subsonic speeds in a constant-area passage under 
the simultaneous influence of friction and heat addition. 

Although the charts are constructed for determining pressure 
drops when heat is added at constant passage-wall temperature, pres- 
sure drops can be determined with good accuracy for other modes of 
heat addition through the use of an effective wall temperature in 
conjunction with these charts. 

The effective wall temperature is given in terms of the passage 
dimensions, gas-flow conditions, and gas-temperature rise across the 
flow passage. The gas-temperature rise is related to the maximum 
passage-wall temperature for the cases of heat addition at constant 
wall temperature, constant rate of heat input along the passage length, 
and sine variation of heat input along the passage length. 


INTRODUCTION 

Many problems arise in aircraft heat-exchanger practice that 
involve the transfer of heat at high flux rates to a compressible 
fluid flowing at high speeds through constant-area flow passages. 

For example, in some recent heat-exchanger designs, the temperature 
differential between heat -exchanger wall and fluid approached values 
of the order of 1200° R and the fluid velocities approached the choke 
condition (Mach number, l). 

Methods have been developed for determining the pressure drop 
experienced in the passage of air through the radiator tubes of 
reciprocating engines (references 1 and 2). These methods contain 
simplifying assumptions that, although adequate over the range of 
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conditions encountered in aircraft radiators, cause appreciable error 
at the high airspeeds and heat-input rates considered herein. 

In reference 3, the basic differential flow equation describing 
the pressure variations of a compressible fluid under the simultaneous 
action of friction and heat addition is numerically integrated for the 
specific case of air heated at constant passage-wall-temperature con- 
ditions. A closed-form solution of the basic differential flow equation 
is obtained in reference 4 for the special case of an exponential vari- 
ation of fluid temperature with distance along the flow passage. In 
reference 5, the basic differential flow equation is formally integrated 
by use of a simplification that is shown to introduce little error in 
the description of the flow process; the integration results are pre- 
sented for air in convenient chart form and are applicable for arbitrary 
heat-input distribution to the fluid along the flow passage. 

Recently, interest has arisen in the use of monatomic gases as 
the working fluid in various aircraft closed cycles. The closed cycles 
involve the use of heat exchangers wherein heat is added to or sub- 
tracted from the monatomic gas. In the calculation of the performance 
of such closed cycles, a convenient means of determining the pressure 
and temperature changes of the monatomic gas flowing through the heat- 
exchanger passages is desirable. 

In this report, which is the result of an investigation made at 
the WACA Lewis laboratory, charts are presented for calculating the 
pressure drop of monatomic gases (ratio of specific heats, 5/3) flowing 
in turbulent motion through smooth constant-area passages wherein heat 
is added to the fluid at constant passage-wall-temperature conditions. 

In the construction of the charts, the recent heat-transfer and fluid 
friction relations obtained in the experimental investigations of 
references 6 to 8 are used. Although the charts are for constant 
passage-wall-temperature conditions, it is shown that through the use 
of an effective passage-wall temperature the charts can be employed 
with good accuracy for other modes of heat input. The effective 
passage-wall temperature is taken as the constant passage-wall tempera- 
ture that results in the same gas -temperature rise as obtained with the 
given mode of heat addition for the same gas-flow conditions and passage 
dimensions . 

Plots are presented that relate the gas -temperature rise to the 
gas-flow conditions, maximum passage-wall temperature, and passage 
dimensions for the following modes of heat addition: 

(a) Constant passage-wall temperature 

(b) Constant rate of heat input along passage length 
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(c) Sine variation of heat input along passage length 

Detailed examples are given that illustrate the use of the charts 
presented herein. 
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SYMBOLS 

The following symbols are used in this report: 
cross-sectional area of flow passage, (sq ft) 


specific heat of gas at constant pressure, (Btu/(lb) (°R) ) 
equivalent diameter of flow passage, which is equal to 4A/s, (ft) 
drag force due to friction, (lb) 


friction factor. 


-d-P 


pV^ 


4dx 


mass conversion factor, 32.2 (lb/slug) 

heat added to gas in distance x per unit cime, (Btu / sec) 

heat added to gas in total tube length L per unit time, 
(Btu) /(sec) 

heat-transfer coefficient between wall and gas, (Btu / 

( sec) ( sq ft) (OR)) 


constant 

thermal conductivity of gas, (Btu/(sec) (sq ft) (°R)/(ft)) 
total tube length, (ft) 

Mach number 

mass flow of gas, (slugs/sec) 
total pressure, (lb/sq ft absolute) 
static pressure, (lb/sq ft absolute) 
gas constant, (ft-lb/(lb)(°R) ) 
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s wetted perimeter, (ft) 

T total temperature of gas, (°R) 

T w temperature of passage walls, (°R) 

T w ^ e effective temperature of passage walls, (°R) 

■*w,max maximum temperature of passage walls, (°R) 

t static temperature of gas, (°R) 

V velocity of gas in flow passage, (ft/sec) 

x distance along flow passage (ft) as measured from station 

where T/T w equals 0.20. This reference station may he a 
hypothetical one located in the forward extension of the 
given flow passage. 


r 

o 

h 

^w,av 


ratio of specific heats 

mass density of fluid, (slugs/cu ft) 

absolute viscosity of fluid, (lb/ft-sec) 

absolute viscosity of fluid evaluated at effective wall tem- 
perature of flow passage, (lb/ft-sec) 


Subscripts: 

en entrance of flow passage 

ex exit of flow passage 

w gas conditions evaluated at passage-wall temperature T w 

f gas conditions evaluated at average of fluid and passage-wall 

T+T w 

temperatures — ^ — 


The following groupings of variables are involved: 

mV + pA 


“p 


total-momentum parameter. 


m 


^gRT - 
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passage-distance parameter 



passage-length parameter 


m 


PA 

\l m 


total -pressure parameter 


PA 

m \fig < T 


static-pressure parameter 


ANALYSIS 

The steady-flow process occurring within a constant- sir ea heat- 
exchanger passage involves the simultaneous action of fluid friction 
and heat transfer. Analysis of this steady-flow process for the case 
of constant passage-wall temperature is presented herein. In addition, 
the applicability of and the method of conveniently using the analysis 
results for modes of heat input other than constant passage-wall tem- 
perature are determined. 

Pressure drop and temperature rise for case of constant wall 
temperature . - One form of the differential momentum equation describ- 
ing the one-dimensional steady-state motion of a compressible fluid in 
a constant-area passage under the combined influence of friction and 
heat transfer is 


d(mV + pA) + dDp = 0 (l) 

From the conservation of energy and mass equations and the perfect gas 
law, it is shown in reference 9 that the total -momentum parameter 

f . P , A (herein denoted NL) is uniquely related to each of the flow 
m^gRT 

parameters M, -, — and — - — for any value of y associ- 

p m^gRT t ^gRT 

ated with the total temperature of the fluid. The equations relating 
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Mp to the aforementioned flow parameters for y equal to a constant 

during the process are presented in convenient form in reference 5 
( appendix B ) . 

It is evident, then, that knowledge of the variations of Mp and 

T during the flow process is sufficient to completely specify the vari- 
ations of all the fluid-flow conditions (p, p, t, and V). 

The variation of Mp during the flow process results from the 
variation of both mV + pA and T. Differentiation of M p with 
respect to these two variables for constant A and m gives 

<M„ - a lrc V +gA) _ 1 (2) 

p -VSt" 2 * T 

The differential drag force dDp is • 

dDp = F | pV 2 A (3) 

or, more conveniently, 

dDp = 2FmV ^5 (4) 

D e 

The variation of T is given by the differential relation equating the 
heat transferred from wall to fluid to the heat absorbed by the fluid 

dH = mgCpdT = hs(T v - T)dx (5) 


The heat-transfer coefficient h and friction factor F are 
evaluated by use of equations (Al) and (A2) of appendix A. These 
equations are based on the recent data of references 6 to 8 obtained 
over a wide range of conditions, including Reynolds numbers up to 
500,000, average tube-wall temperatures up to 2050° R, and heat-flux 
densities up to 150,000 Btu per hour per square foot. A comparison is 
presented in appendix A of pressure drops and temperature rises across 
the test tube as measured in a number of experimental runs and as cal- 
culated by means of equations (Al) and (A2). 
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In view of the analogy between fluid friction and heat transfer 
given by equation (A4) of appendix A for flow through smooth-wall 
passages, equation (5) is rearranged and presented as 


dT = 



( 6 ) 


where, in reducing equation (5) to equation (6), the definition of D e 
and the continuity equation m = pAV are used. 

From equations (l), (2), (4), and (6), the variation of Mp during 
the flow process is expressible as 


dMp = - 


V 


dT 


Y gRT 


1 „ dT 
2^1 


(7) 


CppgV 


F 

2 


|(T W - T) 


The subsonic flow relation between Mp and v/ /^gRT for constant 
X is pre rented in appendix B of reference 5 as 


V 


'll gRT 


r + i 


r4 P 


-{ 


M, 


2 _ 2 (r + 1) 
1 r 


(s) 


Euck n's rule between c ph/k and y, as based on kinetic theory of 
gases, is 


C P^ _ 4 y 

k 9y-5 


( 9 ) 


For V equal to a constant during the flow process, 


T 


w 


dT 


w 


do) 


For monatomic gases , x equals the constant value of 5/3 (hence 
c^ is constant so that in equation (A4) c^ ^/ c v = l) • On the basis 

of equations (A4), (3), (9), and (10), equation (7) reduces to 
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The variation of Mp is uniquely related to the variation of T/T w 
during the flow process by equation (ll). 

In the integration of equation (ll), the lower integration limit for 
T/T„ is arbitrarily taken constant as 0.20, whereas the lower integration 

limit for Mp is assigned about 35 different values that cover the range 
of Mp corresponding to Mach numbers from 0.10 to 0.90. The integration 
is numerically performed to values of Mp corresponding to choke (Mach 
number, l) or to an arbitrarily assigned value of T/T w equal to 0.92, 
whichever is attained first in the integration. 

By equation (ll), M p is shown to be explicitly independent of 
the position variable x; hence each corresponding pair of values of 
T/T w and Mp obtained in an integration performed for a given set of 
lower integration limits can be considered as a possible entrance condi- 
tion of a flow passage. It is evident, then, that although the lower 
integration limit for T/T w is 0.20, the integration results are 

applicable for any values of T/T w at the passage entrance greater 
than 0.20. 

From the integration results and the relations between Mp and 
PA/m^gRT (or pA/m ^gRT) presented in reference 5, the variation 
between PA/m M gRT (or pA/m gP.T) and T/T w is obtained. The 
results are presented in figures 1 and 2 as plots of PA/m ^gRT and 
pA/m ij gRT, respectively, against T/T w . Included in figures 1 and 2 
is a scale giving the Mach number values M corresponding to the 
values of PA/m i\j gRT and pA/m ^gRT. 

The foregoing analysis relates the variation in the fluid-flow 
conditions to the T/T variation obtained during the process. In 

order to obtain the fluid-flow conditions as a function of position 
along a flow passage, it is therefore necessary to know how T/T w 

varies with distance along the flow passage. The following analysis 
provides the required relation between T/T w and x: 
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Substitution of equation (Al) into equation (5) and use of the 
continuity relation ra = pAV and the definition of D e result in 


dT 


0.092 

C P 


kw 


v0.6 




T 


v 0.8 


~ (T V -T, 


^ Y - 2 |2 (12) 

u e 


.mgD, 


(r = 


From equations (9) and (10), for monatomic gases 
5 / 5 ' °P,w/ c P - 11 


ft) 


ftf (■ ■ ft 


= 0.1174 


'M 0 - 2 dx 

mgD ev / 


(13) 


0.1174 


Integration of equation (13) gives 

(Z 2 + 1.6180 Z + l) 0,309 


AmwV * 2 X 

x J 

[rnDj D e 



x o l 


[l-Z)Z 2 - 0.6180 Z + lij 


0.309 


1.9021 tan -1 ( 2Z ~ 0-6180 \ + 1<1756 tan* 1 ( 2Zj_1_ ; _618o' 
V 1.9021 ) [ 1.1756 


(14) 


where Z = (T/T w ) 1//5 . 


Jzq 


The lower integration limits are arbitrarily taken as Xq = 0 
and Z q = (0.20)^/ 3 . These limits are consistent with those taken in 
integration of equation (ll). 

A plot of equation (14) is presented as figure 3 wherein T/T w is 
/Apw \ 0 ’ 2 


the ordinate and 


vmgD e 


is the abscissa. 


Applicability of figures 1 and 2 for other modes of heat input . - 
Although figures 1 and 2 are for the case of constant wall temperature, 
they can be applied to other modes of heat input through the use of an 
effective wall temperature T w ^ e , as shown herein. The effective wall 
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temperature is defined as the constant wall temperature resulting in 
the same temperature rise of the gas for the specified gas-flow con- 
ditions and passage geometry as that obtained with the given mode of 
heat input. Hence, by definition, T w ^ e is the same as T w plotted in 

figure 3. Inasmuch as figure 3 is not in convenient form for determining 
T (or T,, 0 ), a replot of figure 3 is presented as figure 4 where 

W ' W y 6 ' * 


T ex /T en is plotted against 



with T, 


en 


/T„,, 


as the 


curve parameter. 


The following table is presented to show the validity of using 
the effective wall temperature in conjunction with figures 1 and 2 to 
determine pressure drop for modes of heat input other than constant wall 
temperature. Two modes of heat input are considered herein; namely, 
constant rate of heat input, and sine variation of heat input along the 
passage length, both of which represent a radical departure from the 
case of constant passage-wall temperature. The table gives a comparison, 
for these two modes of heat input, of the pressure drop obtained by use 
of figures 1 and 2 in conjunction with the effective wall temperature 
and the pressure drop obtained by integration of the flow differential 
equation for the particular mode of heat input. The comparison is made 
over a range of flow Mach number for T ex /T en equal to 2.0. 


M 

“en 

M 

ex 

Type of heat 
distribution 

ff-) 

venycalc 

(f-) 

V en/ chart 

Difference 

(percent) 

0.207 

0.348 

Constant rate 
of heat input 

0.124 

0.125 

0.8 

.240 

.458 

Constant rate 
of heat input 

.181 

.183 

1.1 

.267 

. 600 

Constant rate 
of heat input 

.244 

.250 

2.5 

.280 

.710 

Constant rate 
of heat input 

.279 

.285 

2.2 

.207 

.348 

Sine variation 
with distance 

.127 

.125 

1.6 

.240 

.458 

Sine variation 
with distance 

.133 

.183 

0 

.267 

. 600 

Sine variation 
with distance 

.250 

.250 

0 

.230 

.702 

Sine variation 
with distance 

.238 

.235 

1.1 
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The agreement indicated in the foregoing table is of the same 
order of magnitude as the accuracy with which figures 1 and 2 can be 
read. It is evident, then, that figures 1 and 2 can be used with 
good accuracy to determine pressure drops in flow passages where 
the heat-input distribution radically departs from the constant 
passage-wall condition by assuming the heat input to the gas to be 
at constant passage-wall temperature with T.. equal to T, 

W - 1 W, 0 


Effective wall temperature for two special modes of heat input . - 
In appendix B, an analysis is presented that relates the temperature 
rise of the gas to the passage dimensions, gas-flow conditions, and 
maximum passage-wall temperature for the cases of constant rate of heat 
input and sine variation of heat input along the passage length. In 
figures 5 and 6, which give the analysis results, T on /T TI . mov is 

Cll' W • IllCvA 


plotted against 


Vw,avY ‘ 2 _L 

k mgD e J D e 


with T ex /T en 


as the curve parameter. 


As demonstrated later in Example III, figures 5 and 6 used in con- 
junction with figure 4 enable convenient determination of T v /T 
for the two modes of heat input treated herein. ' e v ' max 


- 

DESCRIPTION OF CHARTS 


In figures 1 and 2, the ratio of gas total temperature to passage- 
wall temperature T/T v for the case of constant passage-wall tempera- 
ture along the passage length is plotted against the total- and static- 

pressure parameters PA/m^gRT and p*l/m ^gR'T. (in figs. 1 and 2 the 
lines are alternately dashed and solid for convenience in reading.) A 
scale of Mach number to give the values of M corresponding to those 

of PA/m YgRT and pA/m ^gRT is included. The path described by any 
of the curves in figures 1 and 2, going from right to left, describes 
the variation of PA/m^gRT and pA/m /^gRT (or M) with T/T v in a 
constant-area flow passage where T w is constant. Hence, if T/T v 
and M (or PA/m l\J gRT , or pA/m ^g RT) at the flow- passage entrance 
are known, the point on figure 1 or 2 representing the gas conditions 
at the passage entrance is fixed. The curve on which this point lies 
(in fig. 1 or 2) is the line describing the flow variations along the 
passage. If t/t w at the passage exit is known (T ex /T v will be greater 

than T en /T v because of heating) , proceed along the curve to where the 
curve intersects the ordinate value T/T w = T ex /Tw This point 
represents the flow conditions at the passage exit and hence gives the 
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values of M gx and (PA/m^gRT) ex (or (pA/m ^gRT) ex ) . This pro- 
cedure is later illustrated in step-by-step fashion in Example I. 

In figure 3, the ratio of the gas total temperature to the 
passage-wall temperature T/T w is plotted against the distance para- 

\ 0.2 

) for constant passage-wall temperature along the 


meter 


Ap- 


L w 


\fgD e / D e 

passage length. If T/T w at the passage entrance (that is, T en /T w ) 
is known and is equal to or greater than 0.20, figure 3 first gives 

the distance parameter for the passage entrance [that is. 



x is the distance between the passage entrance and a 
en 


Note that 

station in the hypothetical forward extension of the passage character- 
ized by T/T w = 0.20. The distance parameter for the passage exit is 

f Af-L ^ x + L 

then calculated as ( jf- ] e iL where L is the passage length. 


\ingD e y D e 

Figure 3 is then finally used to determine the T/T w corresponding to 
the distance parameter for the passage exit; this value of T/T w is 
the required value of T en /T w - This procedure is later illustrated in 
step-by-step fashion in Example II. 

In figure 4, which is a replot of figure 3 in a form more conven- 
ient for determining T w , is a plot of the ratio of the exit to entrance 
total temper ature of the gas T e x/len against the passage-length para- 
/ A „ \ 0 * 2 

meter | '‘~ H ~ w > av ) _L_ for various ratios of entrance total temperature 




i/T- 


en / w, e' 


As a result of the 


T , figure 4 is applicable for any mode of heat input. 

w j 6 


\mgD e y u e 

of the gas to effective wall temperature T c 
definition of 

The use of figure 4 in conjunction with figures 1 and 2 to determine the 
pressure drop of a monatomic gas flowing in a heated passage for modes 
of heat input other than at constant passage-wall temperature is briefly 
outlined as follows: 

(a) For the given heat-input distribution along the passage 
length, the temperature rise of the gas is determined from the heat- 
balance relation 


nx 


T 


ex 


T 


en 


mgc 


P 


ex 


<3H , 
— dx 
dx 
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where the integration is numerically performed for complicated heat- 
input distributions. 

(b) The value of T w ^ e can be obtained from the known gas-flow 

conditions, passage dimensions, and entrance and exit total tempera- 
tures of the gas by use of figure 4. 

(c) Figures 1 and 2 are then used by treating the problem as one 
for constant wall-temperature conditions with T w = T w g . 

Steps (b) and (c) of the foregoing procedure are illustrated 
later in detail in Example III. 

The relations of the gas total-temperature rise in a flow passage 
to the flow conditions, passage dimensions, and maximum passage-wall 
temperature are given in figures 5 and 6. In figure 5, which is for 
the case of constant rate of heat input along the passage length, the 
ratio of the fluid- entrance temperature to the maximum passage-wall 
temperature T en /T w max is plotted against the passage-length para- 
/ A \0.2 

meter ( av ) J± for various ratios of the exit to entrance total 
\ m S°e / D e 

temperature of the gas T ex /T en . A similar plot for the case of sine 

variation of heat input along the passage length is given in figure 6. 
The use of figures 5 and 6 is described later in detail in Example III. 

For convenience, a plot of the absolute viscosity of helium 
against temperature, as given by reference 10, is presented in figure 7 . 


EXAMPLES ILLUSTRATING USE OF CHARTS 
Example I 

The use of figures 1 and 2 is illustrated for pressure-drop 
determination with heat input at constant passage-wall temperature. 

A monatomic gas is heated during flow through a constant-area 
smooth-wall passage. The passage walls are at a constant temperature 
throughout their length. The flow and heating conditions are as 
follows : 


^en/^w 

T ex/ T w 


( 1 ) 

( 2 ) 

(3) 


0.30 

0.80 

0.20 
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Determine: M 0X and P ex /P en 

(4) The point in figure 1 representing the flow conditions at 
the passage entrance has abscissa (top scale) = 0.20 (item (3)) and 
ordinate = 0.30 (item (l)). This point is designated point a in 
figure 8, which is a small-scale reproduction of figure 1 and is 
presented for illustrative purposes. 

(5) For point a , the lower abscissa scale reads 


P en A 


m 




= 3.974 


en 


(6) From point a , follow the flow-path lines up to the point 
where T/T w = 0.80 (item (2)) in the manner illustrated in figure 8. 
The point where T/T w = 0.80 obtained by this procedure is designated 

point b in figure 8 and represents the flow conditions at the passage 
exit. 


(7) For point b, the lower abscissa scale reads 


P ex A 


1.846 


and the upper abscissa scale reads 


M ex = 0.49 


(8) For T v = a constant, from items (l), (2), (5), and (7)^ 


P A 
ex A 


ex 


■ en 


m 


P A 
r en A 


en 


ex 


w _ 1.346 a h 
T en " 3.974 V‘ 


~w 


0.80 

0.30 


= 0.759 


The procedure for determining the static -pressure ratio P ex /Pen is 
the same as that illustrated for the total-pressure ratio p ex/ p en 
except that figure 2 is used instead of figure 1. 
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Example II 

The use of figure 3 in conjunction with figures 1 and 2 is illus- 
trated for pressure-drop determination with heat input at constant 
passage-wall temperature. 

The flow and heating conditions are as follows: 


(1) T en , °R 500 

(2) T w , °R 1687 

(3) P en , lb/ sq ft 3000 

(4) D e , ft 0.02 

(5) L, ft 1.98 

(6) Mg n 0.20 


(7) Helium gas (R = 386.4) 

Determine: T ex and P ex /P en 

(3) As in Example I, from figure 1 and item (6), 


m Y§RT en 

(9) From items (l), (3), (7), and (8)^ 

ra _ m ^/sE^en P en 
A P en A YgRT en 

1 5000 

3,374 V32.2 X 386.4 X 500 

= 0.303 (slugs/sq ft-sec) 

(10) From figure 7 and items (2) and (l) f 

= 28.1 X io " 6 (lb/ft-sec) 
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(11) From items (4), (5), (9), and (10), 



. 0.2 


D e l 0.303 X 32.2 X o.05 


L / 23.1 X 10" 6 



(12) From items (l) and (2), 


■*-en 500 

T„ " 1667 ~ 0>30 

W 


(13) From figure 3 and item (12), 



which is the distance parameter between a hypothetical station in the 
forward extension of the flow passage, as characterized by T/T w = 0.20, 
and the passage entrance (where T/T w = 0.30 in this case). 

(14) The distance parameter between the hypothetical station 
characterized by T/T w = 0.20 and the passage exit is given simply as 
the sum of items (ll) and (13): 


(16) In order to determine M ex and Pex/Pen> proceed as in 
Example I. Inasmuch as M en , T en /T w , and T ex /T w are the same as in 
Example I, the required values of M ex and P ex /P en 
determined in Example I . 



0.2 


x ex 

= 3.45 + 16.87 = 20.32 


(15) From figure 3 and item (14), 


T 

- L w 


are the same as 
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Example III 


The use of figures 4 and 5 to determine effective wall temperature 
and the application of figures 1 and 2 to the case of constant rate of 
heat input along the flow-passage length are illustrated. 

The heat input to helium flowing in a constant-area smooth-wall 
passage is assumed to he at a constant rate along the passage length. 
The flow and heating conditions are as follows: 


(D T en ,°R 


500 

1900 



16.30 


( 4 ) Men 


0.20 


Determine: T w ^ e , T en /T w>e , T ex /T, 


ex/^w, e* M, 


^ex> p ex/-^en 


(5) From items (l) and ( 2 ) } 


T en _ 500 


= 0.263 


^w,max 1900 


(6) From figure 5 and items (3) and (5), 



(7) From figure 4 and items (3) and (6), 


T 

. en = 0.30 
T 


w,e 


so that from item (l)^ 


= 1667 (°R) 
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(8) From items (6) and (7), 

= o.30 X 2.667 = 0.80 

i w, e 1 w, e i en 

(9) Figures (l) and (2) are used to determine M ex , PexA’en* 
and P ex /p en in the same manner as illustrated in Example I except 
that T/T w ^ e is used in place of T/T w . (Note that in example I, T w 
is a constant, whereas in this example T w is a function of distance 
along the passage. However, T w e is a constant by definition and 
herein corresponds to the constant T w of Example I.) 

For the case of sine variation of heat input, the same procedure 
as outlined in this example is followed except that figure 6 is used 
in place of figure 5. 

Example IV 

The determination of the mass flow of gas through a constant-area 
heated flow passage corresponding to a specified pressure drop is 
illustrated. 


The flow and heating conditions are as follows: 


(1) 

P en' lb / sq ft 


(2) 

P ex' lb / sq ft 


(3) 

^en/^w, e 


(4) 

T /T 

ex' w, e 


(5) 

Ten- 05 



Determine: m/A 


(6) From items (l) to (4), the ratio of the static-pressure 
parameter at the exit and the entrance of the passage must he 


Pex-^ 


m i 




ex 


Pen^ 


m 




en 



en 


w, e _ 3000 
4000 


ex 


w, e 




30 _ 


= 0.459 
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that 


P A 

(7) The problem is to determine the value of - " err such 


p ex A 


m 


/ ygHT en 


(a) 


m 


ViR^ 


ex 


Pen A 


= 0.459 (item (6) ) 


m 




en 


T 

(b) =^L = 0.80 (item (4)) 

W,e 


The procedure used is illustrated in figure 9, which is a small-scale 
reproduction of figure 2 presented for illustrative purposes. 

p A 

(8) Assume — LSI: = 3.50, point a in figure 9. From item 7(a), 

m V^en 

P Q A 

it is necessary that = 0.459 X 3.50 = 1.607. Hence, follow 

m Vg RT ex 

the flow-path line in the manner illustrated up to the point where 
p ex A 

= 1.607, point a in figure 9. 

m VgRT ex 

P A 

(9) Assume en = 4.00 and repeat the procedure outlined 

m^gRT en 

in item (8). Designate the point representing the entrance conditions 
point b and the corresponding point representing the exit conditions 
as point b^ in figure 9. 

(10) Repeat the procedure for an assumed value of 
designating the points in figure 9 points c and c^. 

(11) Draw a curve through the points a}, b^, and 

through the points a, b, and c. The point on curve with 

ordinate value T/T w = 0.80 represents the conditions at the passage 
exit consistent with the given conditions (items (1) through (4)). This 
point is designated point d^ in figure 9. 


^en A 


m 




= 4.50, 


en 


c^ and one 
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(12) The corresponding condition at the passage entrance is 
obtained by following the flow-path lines back from point d^ to 

intersect the curve abc. This intersection is designated point d 

P A 

in figure 9. Point d gives the required value of en satis- 

m ysRT'en 

fying the conditions given in items (7(a)) and (7(b)). 

(13) The lower abscissa value for point d is 


m V gRT en 

(14) Hence from items (l), (5), and (13) for R = 386.4 (helium) 

m _ m V^en Pen 
A Pen A ^iRT^ 


m = _1 4000 

A 4.282 y\/32.2 X 386.4 X 500 


0.374 ( slugs/sec-sq ft) 


Lewis Flight Propulsion Laboratory, 

National Advisory Committee for Aeronautics, 
Cleveland, Ohio, November 3, 1950. 
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APPENDIX A 

HEAT-TRANSFER AND FLUID -FRICTION RELATIONS 

Some recent data available on heat-transfer coefficient and 
friction factor for fluid flow through smooth passages are given by 
the extensive experimental investigations of references 6 and 7. 

These data were obtained with air flowing through an electrically 
heated Inconel tube with an inside diameter of 0.4 inch and a length 
of 24 inches for various types of tube-entrance section. The range of 
conditions investigated included Reynolds numbers up to 500,000, average 
tube-wall temperatures up to 2050° R, and heat flux densities up to 
150,000 Btu per hour per square foot. Based on the correlation results 
of references 6 to 8 obtained for the entire range of surface tempera- 
tures and for Reynolds numbers above 10,000, the following equation 
for local heat-transfer coefficient h is used in the analysis of this 
report: 



A comparison, for a number of experimental runs, between experimental 
and calculated temperature-rise values (using equation (Al)), gave 
agreement within 10 percent. 

Several equations for correlating friction factor, none of which 
is completely satisfactory, were used in references 7 and 8; neverthe- 
less, improved correlation over that given by the conventional correl- 
ation method was obtained. The two more satisfactory correlation 
equations are given by 


F = 



(A2) 



(A3) 


The validity of using equation (a 2) or (A3) for calculating local 
friction factors is measured by how accurately pressure-drop predictions 
can be made through the use of these equations. Comparison is made herein 
between calculated and measured pressure variations along the test passage 
of references 7 and 8 for a number of experimental rims conducted at high 
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heat inputs to the fluid at both low- and high-flow Mach numbers. 

The calculated pressure variations were determined using both equa- 
tions (A2) and (A3) in order to determine which of these two equations 
results in best agreement with the experimental data. The test condi- 
tions of the experimental runs and the comparison results are tabulated 
as follows: 


M 

en 

M 

ex 

T 

w,av 

T 

en 

T 

1 ex 

Ap/Pen 

Experi- 

mental 

Calculated 
(equation (A2)) 

Calculated 
(equation (A3)) 

0.410 

0.722 

1072 

554 

753 

0.361 

0.348 

0.418 

.400 

.770 

1273 

530 

798 

.346 

.326 

.392 

.358 

.732 

1677 

525 

917 

.381 

.350 

.433 

.411 

.836 

1196 

543 

781 

.352 

.326 

.392 

.172 

.234 

1246 

552 

904 

. 063 

. 058 

.063 

.374 

.839 

1679 

528 

903 

.333 

.352 

.425 

.353 

.632 

1272 

529 

312 

.316 

.269 

.317 


Inasmuch as in the experiments the heat was added to the fluid at 
essentially constant rate along the passage, the fluid temperature was 
taken to vary linearly across the passage. In order to reduce entrance 
effects, the fluid conditions to start the calculations were taken as 
those existing in the test passage at a distance of 4 inches (10 diam.) 
downstream of the entrance. 


The foregoing table indicates that equation (A2) gives better 
agreement with the experimental data at the high Mach numbers; at the 
low Mach numbers, equation (A3) gives better agreement. 

For the analysis presented in this report, equation (A2) is chosen 
to represent the friction process. 

From equations (Al) and (A2), 


h 

SpVg 

F 

2 



c p,wt-*w 



(A4) 
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APPENDIX B 

RELATION BET 1 , -TEEN GAS -TEMPERATURE RISE AND MAXIMUM WALL TEMPERATURE 
FOR CONSTANT HEAT FLUX ALONG PASSAGE AND SINE HEAT-INPUT 
VARIATION ALONG PASSAGE 

Relations are derived herein that will permit determination of 
the temperature rise of the gas in terms of the passage dimensions, 
gas-flow conditions, and maximum passage-wall temperature for the 
cases of constant heat flux along the passage and sine heat-input vari- 
ation along the passage. 


Constant heat flux along passage . - The case of constant heat flux 
is characterized by 


dH 

dx 


H+ 


(Bl) 


From equations (5) and (Bl), the continuity equation m = pAV, and the 
definition of D e , 




4h / ^w 


(L/D e )mgc T pVgc^ V T 


- 1 


(B2) 


Subs Ltution of equation (Al) in equation (B2) results in 


1 - / k „ 

7T-7TT rn = 0.032 

mgc p T l° p ,v*V 



(B3) 


wherein, f >r monatomic gases, c p,w/ c p = 

The t( tal amount of heat added to a fluid in the over-all passage 
length L can be given as 


H t = mgc (T e -T en ) 


(B4) 
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At any flow cross section of the passage, equation (B3) is appli- 
cable. When applied to the exit of the passage, T in equation (B3) 
is replaced by T gx and subscript w is replaced by subscript w,ex. 
Hence, from equations (9), (B3), and (B4) for y = 5/3, 



For the case of constant heat flux along the passage, the exit passage- 
wall temperature T w ^ ex is the maximum wall temperature T w , max • A 

plot of equation (B5) is presented as figure 5. Because 

varies only slightly with temperature, is replaced by u in 

figure 5 in order to permit convenient use of figure 5 in conjunction 
with figure 4 for determining effective wall temperature. 

Sine heat-input variation along passage . - For the case of sine 
heat-input variation along the passage. 


Combination of the relation dH = mgc dT (equation (5)) and equation (B6) 

Jr 

and integration result in 


dH = K sin — dx 
L 


(BS) 



(B7) 


In equation (B7), when x = L, T = T gx so that 



(B8) 
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Combination of the relation dH = hs(T w -T)dx (equation (5)) and 
equation (B6) and substitution for h from equation (Al) give 



(B9) 


where 


r = 0.023 Cp^O-2 



(Bio) 


The condition of maximum wall temperature is given by dT w /dx = 0 

Differentiating equation (B9) with respect to x and setting 
dT w /dx = 0 give 




cos 


rtx 


L 

YL dT 
Krt dx 


0 


(Bll) 


Substitution of equation (B9) into equation (Bll) and use of the 
relation ^ . sin (given by equations (5) and (B6)) result 

in P 



„ o Jtx 
0.8 sin -j- 


YT 

K 


. 7TX 

cot — 
YL 

mgc p * 


(B12) 


From the definition of Y, let 

_ YL _ 0.09 2 f \°’ 6 A J v,av\°* 2 L 
rtmgC p" * ^ c p,w^ [ m ^eJ D e 


(B13) 


Note that 


YT = 
K 


range T 


KL 


E— A, which; from equations (B7) and (BQ), gives 


K 

K 




(B14) 


26 


NACA TN 2328 


Hence, equation (B12) may be written as 



(B15) 


where $ is a function of x/L and T ex /T en as defined by equa- 
tion (B15) . 


Substitution of equations (B14) and (B15) in equation (B9) and 
solution for A give 

'| 0.8 


A = 


sm 


nx 


1 + 




2 

— 


^ex n ( 

n *x\ 


T ” 1 + l 

en \ 

1 - COS -) 

cp 






qp 


(B16) 


In equations (B15) and (B16), A and qp are given as functions of the 
basic variables x/L and T ex /T en . From equation (B15) , the corres- 
ponding value of T/T , which is given by dT w /dx = 0 and hence is 

actually T/T tt is determined. From equations (B7) and (B8), 

1 w y max 

T^„/T, r is then obtained. Note that for monatomic gases (y = 5/3 

6 n/ w y max 

and c-rWk = 2/3 by equation (9)), A = 0.03736 [ av ) — by 

p V mgD p / D p 

equation (B13). A plot of T e n/ T w,max agai- 11 ^ f —~ — ) — por 

^ mgD e J D e 

various values of T ex /T en is presented as figure 6. 
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Figure 1. - Variation of total-pressure parameter with ratio of fluid total temperature to passage- 
wall temperature. Turbulent flow of monatomic gases (7 = 5/3) through smooth passages with 
constant wall temperature. (A 22- by Q^-in. print of this fig. is attached.) 
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Figure 3. - Variation of ratio of fluid to wall temperature with distance 
parameter for monatomic gases (7 = 5/3) and constant wall-temperature 
conditions. (A 19- by 17 -in. print of this fig. is attached.) 
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(a) Ratio of exit to entrance fluid temperature, 1.0 to 2.02. 

Figure 4. - Chart for determination of effective constant wall temperature that 
results in any given fluid-temperature rise within flow passage. Turbulent 
flow of monatomic gases (y = 5/3). A 15- by 21- in. print of this fig. is 

attach'd ) 
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Figure 4. - Concluded. Chart for determination of effective 
constant wall temperature that results in any given fluid- 
temperature rise within flow passage. Turbulent flow of 
monatomic gases (y = 5/3). a 11- by 15- in. print of this 
fig. is attached.) 
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’ Var * atl °" of ratl ° of fluid-entrance temperature to maximum passage-wall temperature 
with passage-length parameter and fluid-temperature ratio. Turbulent flow of monatomic gases 
[y - 5/3} for conditions of constant rate of heat input. (A 22- by 16- in. print of this fig 
is attached. ) 
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Figure 6. - Variation of ratio of fluid-entrance temperature to maximum passage-wall temperature with 
oas sage -length parameter and fluid-temperature ratio for sine variation of heat input to fluid. 
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Figure 7. - Variation of absolute viscosity of helium with temperature 


Choke line, PA/m^/gRT = 1.377 
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